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Abstract Oryza officinalis complex includes five diploid species involving the B, C and E 
genomes, and is distributed in Asia, Africa and Oceania. These species are very important because 
O. australiensis is the only species with E genome, while O. punctata is the single species with B 
genome. Although there are three species with the C genome (O. eichingeri, O. officinalis and 
O. rhizomatis), these species might have played different roles in speciation of polylploids due to 
their disjunction in Asia and Africa. The phylogenetic relationships among these diploids are still 
unclear, and polyploidy speciation involving these diploids has long been in debate. In this study, 
we utilized the sequence data of multiple genes, including the chloroplast matK, internal tran- 
scribed spacer (ITS) fragment of nuclear ribosomal DNA, and two nuclear genes ( Adhl and 
Adh2), to infer the phylogenetic relationships among the diploids and the related genomes. The 
separate and combined analyses were implemented by both maximum parsimony and maximum likeli- 
hood methods. The results indicate that (1) the C genome showed a closer affinity to the B genome 
than to the E genome; (2) of the three species with the C genome, O. officinalis and O. rhizoma- 
tis showed a closer relationship to each other than to O. eichingeri; (3) notably, in O. eichingeri , 
there exists higher divergence between the samples from Africa and those from Sn Lanka although 
they formed one monophyletic clade with lower bootstrap support. 

Key words phylogenetic relationships, Oryza officinalis complex, diploid, multiple gene 
sequence . 

摘要 ” 稻 属 药 称 复合 体 Oryza officinalis complex 中 有 5 个 二 倍 体 物 种 ,涉及 B.C 和 了 3 个 染色 体 组 ,分 布 
在 亚洲 ,非洲 和 大 洋 洲 。 其 中 ，0O .australiensis 是 舟 属 中 惟一 含 了 染色体 组 的 物种 ; 0. punctata 则 是 惟 
一 含 B 染 色 体 组 的 物种 。 虽 然 0 officinalis, O. rhizomatis 和 O. eichingeri 都 含 C 染色 体 组 ,但 由 于 它们 
间断 分 布 在 亚洲 和 非洲 ,可 能 在 稻 属 异 源 多 倍 体 成 种 中 发 挥 了 不 同 作用 。 染 色 体 组 B`C 和 下 之 间 ,以 及 
上 述 二 倍 体 物种 之 间 的 系统 发 育 关系 及 其 在 稻 属 多 倍 体 物种 形成 中 的 作用 都 是 稻 属 系统 和 进化 研究 中 
长 期 争论 的 问题 。 本 文 利用 多 基因 片段 ,包括 叶绿体 matK 基因 、nrDNA ITS 片段 、 低 拷贝 核 基 因 Adhl 和 
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Adh2, 探 讨 了 这 些 二 倍 体 种 以 及 它们 所 代表 的 染色 体 组 之 间 的 系统 发 育 关系 。 独 立 和 合并 的 分 析 显 示 
了 一 致 的 系统 发 育 关 系 , 即 C 染色 体 组 和 B 染 色 体 组 的 亲缘 关系 要 比 它们 和 下 染色 体 组 的 更 近 。 在 3 
个 含 C 染 色 体 组 的 二 倍 体 中 , O. officinalis WO. rhizomatis 表现 出 较 近 的 亲缘 关系 。 值 得 注意 的 是 ,在 
O. eichingeri 种 内 ,尽管 基于 多 基因 的 数据 支持 来 自 斯 里 兰 卡 的 样品 和 来 自 非 洲 的 样品 为 单 系 分 支 ,但 
是 较 低 的 支持 率 表明 这 两 个 地 区 的 样品 间 存 在 着 较 大 的 遗传 分 化 ,值得 进一步 利用 形态 .遗传 和 地 理学 
证 据 从 居 群 水 平 上 进行 深入 的 研究 。 

关键 词 ”系统 发 育 关系 ; 药 稻 复 合体 ; 二 倍 体 ; 多 基因 序列 


药 称 复合 体 ( Oryza officinalis complex) 是 稳 属 中 最 大 的 复合 体 ,包括 9 个 种 , 含 5 种 染色 
体 组 类 型 (B,C,BC,CD,E)(Vaughan，1989，1994), 广 泛 分 布 于 亚洲 非洲 、 大 洋 洲 和 拉丁 
美洲 。 在 该 复合 体 中 ,5 个 二 倍 体 种 在 稻 属 中 具有 非常 重要 的 地 位 。 其 中 , O. australiensis 
Domin .是 稻 属 中 惟一 含 下 染色 体 组 的 物种 ; O. punctata Kotechy ex Steud. 则 是 惟一 含 B Ye 
色 体 组 的 二 倍 体 和 种。 虽然 O. officinalis Wall. ex Watt., O. rhizomatis Vaughan 和 O. 
eichingeri A. Peter 均 含 C 染色 体 组 ,但 它们 间断 分 布 在 亚洲 和 非洲 ,在 稻 属 异 源 多 倍 体 
BBCC 和 CCDD 的 物种 形成 中 具有 不 同 的 作用 。 但 是 ,这 些 二 倍 体 物 种 之 间 的 系统 发 育 关 
系 以 及 一 些 物种 的 分 类 地 位 一 直 存 在 争议 (Bor，1960; Sharma & Shastry, 1965; Tateoka, 
1962; Vaughan，1990)。 例 如 ,根据 Vaughan (1990) 的 研究 ,斯里兰卡 只 有 2 个 含 C 染色 体 
组 的 二 倍 体 种 , 即 O. rhizomatis FO. eichingeri ,前 者 是 斯 里 兰 卡特 有 种 ,后 者 则 间断 分 布 
于 亚洲 和 非洲 ,在 亚洲 也 只 局 限于 斯 里 兰 卡 (Vaughan，1994)。O.， eichingeri 曾 被 定名 为 
O. collina (Sharma & Shastry, 1965) ,后 Vaughan(1990) 将 其 并 入 O. eichingeri ,但 对 这 些 物 
种 的 分 类 仍 缺 乏 统一 的 认识 。Wang 等 (1992) 利 用 RFLP 研究 稻 属 物种 的 多 态 性 和 亲缘 关 
系 时 ,将 编号 为 103410 和 103421 的 2 份 材料 定名 为 0. collina ,而 Vaughan 1990 年 发 表 新 
种 O. rhizomatis 时 已 将 这 2 份 材料 处 理 为 0. rhizomatis。 由 于 不 同学 者 在 以 往 的 分 类 处 
理 时 基本 依据 形态 特征 和 个 人 经 验 ,因此 迫切 需要 利用 新 的 技术 ,特别 是 分 子 技术 探讨 这 
些 二 倍 体 的 系统 发 育 关 系 ,为 物种 界定 提供 进一步 的 证 据 。 同 时 ,探讨 药 稻 复合 体 二 倍 体 
种 间 的 系统 发 育 关 系 对 于 深入 理解 整个 称 属 的 起 源 和 分 化 以 及 物种 形成 方式 具有 重要 意 
Xo 

叶绿体 DNA(cpDNA) 和 细胞 核 核 糖 体 DNA (nrDNA) 是 系统 学 研究 中 最 常用 的 分 子 标 
记 。 尽 管 大 多 数 叶 绿 体 基因 进化 速率 较 慢 ,从 而 主要 用 于 探讨 较 高 级 分 类 等 级 的 系统 发 
育 关系 ,但 一 些 进 化 稍 快 的 叶绿体 基因 (如 matK) 或 基因 间隔 区 (如 wmL-F) 等 ,被 有 效 地 
用 于 较 低 等 级 的 系统 发 育 研 究 (Johnson & Soltis, 1994; Sang et al., 1997; Bakker et al. , 
1999; Kajita et al., 1998; Ge et al., 1999; Wang et al., 1999), HH nrDNA 的 转录 间隔 区 
(ITS) 是 经 历 了 致 同 进 化 (concerted evolution ) 的 多 拷贝 片段 ,进化 速率 快 且 易于 PCR 扩 增 
和 测序 ,从 而 被 广泛 地 用 于 低 等 级 的 系统 发 育 研 究 中 ,对 揭示 近 缘 属 和 属 内 种 间 的 系统 发 
育 关系 提供 了 重要 的 依据 (Baldwin, 1992; Sang et al., 1997; Baldwin et al., 1995; Soltis & 
Soltis，1998)。 近 年 来 ,一 些 低 拷贝 核 基 因由 于 包含 多 个 内 含 子 和 外 显 子 ,含有 丰富 的 系 
统 发 育 信息 ,被 成 功 地 应 用 于 各 个 不 同等 级 的 系统 发 育 重建 ,并 显示 了 巨大 的 潜力 (Soltis 
& Soltis, 1998; Ge et al., 1999; Sang，2002)。 但 是 , 单 基因 片段 提供 的 系统 发 育 信息 有 
限 ,而 且 采 用 单 基因 片段 构建 的 系统 树 可 能 反映 的 仅仅 是 该 片段 的 基因 树 ,在 许多 情况 下 
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基因 树 并 不 等 同 于 物种 树 (Soltis et al., 1996; Wendel & Doyle, 1998; Sang,2002)。 因 此 ， 

当前 分 子 系统 发 育 重建 的 一 个 主要 趋势 是 通过 多 个 独立 片段 构建 系统 发 育 树 并 进行 比较 
分 析 , 最 后 推断 物种 的 系统 发 育 关系 (Hillis,，1995; Wendel & Doyle, 1998; Soltis & Soltis, 
1998; Sang,，2002)。 本 研究 选取 了 叶绿体 matK 基因 nrDNA 的 ITS 片段 和 低 拷贝 核 基 因 
乙醇 脱氧 酶 基因 Adhl 和 4dh2。 通 过 这 些 DNA 片段 的 独立 和 合并 分 析 ,探讨 稳 属 药 称 复 
合体 中 二 倍 体 物种 间 的 系统 发 育 关 系 ,为 该 复合 体 以 及 稳 属 的 分 类 提供 新 的 资料 ,同时 为 
进一步 阐明 B.C 和 下 染色体 组 之 间 的 关系 以 及 揭示 该 复合 体 中 多 倍 体 物种 形成 方式 和 
机 制 莫 定 基础 。 


1 材料 和 方法 


1.1 材料 

本 研究 所 用 的 材料 代表 药 稻 复合 体 的 5 个 二 倍 体 物种 。 每 个 二 倍 体 种 原则 上 选取 2 
份 样品 。 由 于 Oryza eichingeri 间断 分 布 于 亚洲 和 非洲 , 故 选取 了 4 份 材料 ,2 份 来 自 非洲 ， 
2 份 来 自 斯 里 兰 卡 ; 仅 分 布 于 澳大利亚 的 O，australiensis 除 在 ITS 片段 分 析 时 选用 了 2 份 
样品 ,其 他 基因 片段 均 选 用 1 份 样品 。 男 外 ,1 份 来 自 中 国 海南 的 O. granulata Nees & 
Am. ex Watt 样品 被 选 作 外 类 群 ,因为 该 种 被 证 明 是 稳 属 最 基部 的 谱系 (Wang et al.， 
1992; Ge et al.，1999)。 上 述 12 份 样品 的 详细 资料 见 表 1。 
表 1 植物 样品 材料 


Table 1 Plant materials used in the present study 





el, weep EE ers 
样品 染色 体 组 。 IRRI 编 号 原 产 地 序列 接收 号 GenBank numbers 
Sample!” Genome IRRP? accession Origin matK ITS Adh1 Adh2 
Oryza officinalis -1 cc 105085 菲律宾 Philippines AF148658 AF479067 AF148579 AF148613 
O. officinalis-2 CC 106159 巴布亚 新 几内亚 - AF479063 AY169473 AY169483 
Papua New Guinea 
O. eichingeri-sl cc 81803 斯 里 兰 卡 Sri Lanka AY176644 AF479068 AY169475 AY169487 
O. eichingeri-s2 CC 105407 斯 里 兰 卡 Sri Lanka - AY181994 AY169476 AY169486 
O. eichingeri-ul CC 101422 乌干达 Uganda AY318858 AF479069 AY169477 AY169485 
O. eichingeri-u2 cc 105159 乌干达 Uganda - AY181995 AY169478 AY169488 
O. rhizomatis-1 CC 105448 斯 里 兰 卡 Sri Lanka AF148660 AF479066 AF148580 AF148614 
O. rhizomatis-2 CC 103410 斯 里 兰 卡 Sri Lanka - AF479065 AY169474 AY169484 
O. punctata-1 BB 104071 乍得 Chad AF148661 AF479070 AF148577 AF148611 
O. punctata-2 BB 105607 乍得 Chad - AF479071 - - 
O. australiensis-1 EE 105263 澳大利亚 Australia AF148667 AF520777 AF148589 AF148623 
O. australiensis-2 EE 101144 澳大利亚 Australia - AF520778 - - 
O. granulata 940208 中 国 海南 Hainan, China AF148674 AF520779 AF148597 AF148631 


1) eR s 代 表 斯 里 兰 卡 ,代表 乌干达 。2) IRR, 国际 水 称 研 究 所 。 
1) Numbers after species name indicate different accessions; s indicates Sri Lanka; u indicates Uganda. 2) IRRI, International Rice 
Research Institute. 


除外 类 群 为 野外 采集 ,其 余 的 材料 均 来 自 国际 水 称 研 究 所 (IRRI) 的 种 子 资源 库 。 野 
生 稻 种 子 首先 在 50 - 55 C 处 理 5 d 以 打破 休 阴 ,然后 将 去 壳 的 种 子 悬 浮 在 灭 菌 液 中 充分 
冲洗 ,在 温水 (30 - 35 °C) PARA. ÆR 2 2 周 的 幼苗 被 移植 至 中 国 科学 院 植物 研究 所 ( 北 
京 ) 的 温室 中 栽培 。 用 2 个 月 幼苗 的 叶片 提取 DNA。 取 样 植株 均 保留 至 成 熟 , 以 观察 其 形 
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态 特 征 。 
1.2 DNA 的 提取 、PCR 扩 增 和 DNA 测序 

从 个 体 新 鲜 叶 片 中 提取 总 DNA 均 使 用 改良 的 CTAB 法 (Doyle & Doyle, 1987). RA 
酶 链 式 反应 (PCR) 在 PTC-200(PE) 热 循环 仪 上 完成 。 扩 增 叶 绿 体 matK 基因 的 程序 包括 
94 CHAHE 1.5 min;35 个 循环 的 94 CASPE 1.5 min,55 CIRK 1.5 min 和 72 延伸 1.5 
min; 最 后 72 所 延伸 5 ming ITS 的 扩 增 反应 中 加 入 8% 的 二 甲 亚 砚 (DMSO) 以 增加 引物 特 
异性 (Baldwin et al., 1995) , 扩 增 程序 包括 70 CHA 4 min; 2 个 循环 的 94 CAEPE 1 min, 
52 CiB-K 20 s 和 72 % 延 伸 1.5 min; 接 35 个 循环 的 94 CARTE 20 s,52 CIBK 20 s 和 72 
延伸 1.5 min; 最 后 72 延伸 10 mino Adhl 或 Adh2 基因 的 扩 增 程序 包括 70 °C HA 4 
min;2 个 循环 的 94 变性 1 min,56 TC 退火 30 s 和 72 CHE 1.5 min; 接 35 个 循环 的 94 
GC 变性 20 s,55 CIRK 20 s 和 72 CHEM 1.5 min; 最 后 72 CEH 10 min。 扩 增 及 测序 所 
用 的 引物 见 表 2。 扩 增产 物 利 用 Pharmacia 纯化 试剂 盒 纯化 后 在 ABI377 自动 测序 仪 上 进 
行 测序 。 为 保证 序列 的 精确 性 ,所 有 序列 均 对 双 链 进行 测序 。 


表 2 用 于 扩 增 和 测序 的 引物 序列 
Table 2 The primer sequences for PCR amplification and sequencing 


引物 序列 参考 文献 
Primer Sequence Reference 
mat K Ge et al., 1999 
matKF 5'-TAATTAAGAGGATTCACCAG-3’ 
matKF2 5'-ATTGCCTTTCCTTGATATCG-3’ 
matKR2 5’-ACTACTCGAATTGGAATAG-3' 
matKR1 5'-ATGCAACACCCTGTTCTGAC-3' 
ITS White et al. , 1990 
ITS1 5'-AGAAGTCGTAACAAGGTTTCCGTAGG-3' 
ITS4 5'-TCCTCCGCTTATTGATATGC-3’ 
Adh1 Ge et al. , 1999 
AdhF 1 5’-CACACCGACGTCTACTTCTG-3’ 
AdhF2 5'- AGAGTGTTGGAGAGGGTGTGAC-3’ 
AdhI R2 5’-ACTCACAGCAAGGCCTACAGC-3’ 
AdhlbR 5'-TCAGCAAGTACCTAAATTATC-3' 
Adh2 Ge et al. , 1999 
AdhF1 5’-CACACCGACGTCTACTTCTG-3’ 
AdhF2 5’-AGAGIGITGGAGAGGGTGTGAC-3’ 
Adh2R2 5’- ACAGCAAGGCCAACAGCTCC-3’ 
Adh2RR 5’-CCACCGITGGTCATCTCAAT-3’ 
1.3 数据 分 析 


序列 通过 CLUSTAL W (Thompson et al. ，1994) 进行 对 位 排列 (alignment) 后 手工 校正 ; 
系统 发 育 分 析 采 用 PAUP 4.0 版 (Swofford, 1998) 中 的 简约 法 和 似 然 法 进行 。 最 大 简约 性 
(MP) 分 析 选 择 分 支 限 界 (branch and bound search) 法 进行 ,对 所 有 核 苷 酸 位 置 的 碱 基 性 状 
同等 加 权 (Fitch，1971)。 对 于 由 于 序列 长 度 多 态 性 所 造成 的 空位 (gap) ,在 数据 统计 中 不 
参与 运算 , 即 处 理 为 缺失 (missing) 状 态 。 最 大 似 然 性 (ML) 分 析 依 据 HKY 序列 进化 模式 
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(Hasegawa et al., 1985) ,选择 分 支 限界 法 进行 ML 搜寻 1000 次 。 碱 基 频 率 设 定 为 数据 佑 
$8. (empirical) ,转换 和 其 换 比 设 定 为 2。 拓扑 结构 的 可 靠 性 用 1000 次 重复 的 自 展 (beot- 
strap) 检 验 来 完成 (Felsenstein，1985)。 为 了 检测 来 自 不 同 DNA 片段 所 得 到 的 系统 发 育 树 
是 否 存 在 冲突 ,利用 PAUP 4.0 版 程序 中 的 partition-homogeneity test (PHT) 对 叶绿体 matK 
基因 、 核 ITS H Bt Adh 1 和 Adh2 基因 系统 树 之 间 的 拓扑 一 致 性 进行 评价 (Johnson & Soltis, 
1998) ,检测 进行 1000 次 重复 。 当 来 自 不 同 DNA 片段 所 得 到 的 系统 发 育 树 间 不 存在 显著 
冲突 (P>0.05) 时 ,不 同 片 段 的 数据 将 进一步 合并 进行 系统 发 育 分 析 。 


2 结果 


2.1 序列 特点 

扩 增 的 叶绿体 matK 基因 片段 长 1487 - 1552 bp。 对 位 排列 后 的 长 度 为 1552 bp, 包 括 
34 个 变异 位 点 ,其 中 15 个 含 系统 发 育 信息 。TTS 片段 的 长 为 578 - 610 bp, 包括 ITS1(188 
— 220 bp) ,5.8S(164 bp) 和 ITS2(228 - 230 bp)3 个 区 段 ,对 位 排列 后 的 长 度 为 668 bp, 含 56 
个 变异 位 点 ,其 中 33 个 含 系统 发 育 信息 。 尽 管 TTS1 的 长 度 短 于 ITS2, 但 含 系统 发 育 信息 
的 位 点 数 (16 个 ) 与 ITS2(16 个 ) 相 同 ,5.8S 基因 序列 高 度 保守 ,在 5 个 二 倍 体 物种 中 只 发 
现 2 个 变异 位 点 ,其 中 1 个 含 系统 发 育 信息 。 扩 增 的 Adhl 和 Adh2 基因 片段 包括 5 个 内 
含 子 和 5 个 外 显 子 。 扩 增 的 4dh1 基因 片段 长 1770 - 1829 bp, 对 位 排列 后 的 长 度 为 1940 
bp, 包 含 264 个 变异 位 点 ,其 中 133 个 含 系统 发 育 信 息 ; Adh2 片段 长 1645 - 1668 bp, 对 位 
排列 后 的 长 度 为 1752 bp, 包 含 236 个 变异 位 点 ,其 中 137 个 含有 系统 发 育 信息 。4 个 片段 
中 , 核 基 因 Adhl 和 Adh2 提供 的 系统 发 育 信 息 最 多 (信息 位 点 分 别 占 7% 和 6%);ITS 次 之 
(5%) ;叶绿体 matK 片段 提供 的 信息 最 少 (1% )。 另 外 ,TTS 片段 的 GC 含量 (72% ) 要 明显 
高 于 其 他 片段 (35% - 43% )( 表 3)。 


表 3 各 DNA 片段 的 序列 参数 
Table 3 Sequence characteristics for each DNA fragment 


序列 排列 长 度 。 变异 位 点 。 “信息 位 点 信息 位 点 (%) PK 





apie at Aligned Variable Informative Informative Tree CI RI Gs 
length (bp) site site site (% ) length 

mat K 1552 34 15 1 425 0.93 0.90 35 

ITS 668 56 33 5 57 0.97 0.91 72 

Adh1 1940 264 133 7 110 0.85 0.84 4l 

Adh2 1752 236 137 6 392 0.93 0.90 43 

mat K/ITS/ Adh 1 4160 354 181 4 590 0.92 0.88 44 


2.2 系统 发 育 分 析 

对 各 个 DNA 片段 分 别 采 用 最 大 简约 法 (MP) 和 最 大 似 然 法 (ML) 进 行 系统 发 育 分 析 。 
4 个 基因 片段 各 自得 到 最 简约 树 的 数目 分 别 为 1 个 (matK)、2 个 (ITS)、2 个 (44dh1) 和 3 个 
( 4Adh2) ,其 严格 一 致 树 的 树 长 .一致 性 指数 ( C7) 和 保持 性 指数 (RD) 列 于 表 3。ML 得 到 的 
系统 发 育 树 在 分 支 结 构 上 与 MP 的 结果 一 致 , 仅 在 自 展 支持 率 上 略 有 不 同 , 故 在 此 只 展示 
ML 的 结果 。 图 1 - 3 为 根据 ML 分 别 得 到 的 4 个 DNA 片段 的 单 基因 系统 发 育 树 。 由 图 可 
见 ,4 个 ML 树 都 是 按照 染色 体 组 CB 和 了 分 成 3 个 主要 分 支 , 含 C 染色 体 组 物种 构成 的 
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分 支 和 含 B 染色 体 组 物种 构成 的 分 支 形 成 姐妹 群 关系 ,然后 再 与 含 下 染色 体 组 物种 的 分 
文 构成 姐妹 群 。 这 4 个 基因 树 的 主要 不 同 在 于 C 染色 体 组 分 支 内 部 。 在 muk 基因 树 
EL, 0. officinalis 的 不 同样 品 形成 一 单独 小 分 支 (支持 率 81% ) ,其 余 两 个 二 倍 体 以 多 歧 分 
文 (polytomy) 的 形式 出 现 (图 1)。 然 而 在 ITS 的 ML 树 上 , O. officinalis O. rhizomatis 都 
各 自 构 成 单独 的 小 分 支 ,而 O. eichingeri 的 斯 里 兰 卡 样品 和 乌干达 样品 构成 了 另外 两 个 
小 分 支 , 这 两 个 小 分 支 构成 的 单 系 未 得 到 显著 的 支持 ( < 50% )( 图 2)。 在 hdh1 的 基因 树 
上 ,0O. officinalis 表现 出 和 O. rhizomatis 更 近 的 关系 (支持 率 81% ) AA O. eichingeri 的 2 
个 地 区 样品 也 没有 构成 单 系 (图 3); 4dh2 的 ML 树 的 基本 结构 和 Adh1 的 类 似 , 主要 冲突 
在 于 O. rhizomatis 的 两 个 样品 没有 构成 独立 的 小 分 支 ,其 中 一 个 样品 先 和 O. officinalis 
小 分 支 聚合 ,再 和 另 一 个 样品 形成 姐妹 群 (图 3)。 另 外 ,这 4 个 基因 树 在 分 辨 率 上 也 存在 
不 同 , 其 中 mat Kk 基因 树 的 分 辩 率 最 小 ,ITS 的 次 之 ,4 基因 的 最 大 ,体现 了 各 DNA 片段 
所 含 信 息 量 的 差异 。 
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图 1 叶绿体 mark 基因 的 最 大 似 然 树 (种 名 后 的 数字 
代表 不 同样 品 ,s 代表 斯 里 兰 卡 ,u 代表 乌干达 ;分 支 上 
eas 50% 的 自 展 支持 率 ,字母 代表 染色 体 组 类 


Fig. 1. The maximum likelihood tree based on matK. Num- 
bers after species name indicate different accessions; s indi- 
cates Sri Lanka; u indicates Uganda; numbers above branches 
indicate bootstrap values above 50% ; letters in capital repre- 
sent the genome types. 














2.3 不 同 DNA 片段 的 合并 分 析 
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图 2 ITS 的 最 大 似 然 树 (种 名 后 的 数字 代表 不 同样 品 ,s 
代表 斯 里 兰 卡 ,u 代表 乌干达 ;分支 上 的 数字 为 大 于 
50% 的 自 展 支 持 率 ,字母 代表 染色 体 组 类 型 ) 

Fig. 2. The maximum likelihood tree based on ITS. Num- 
bers after species name indicate different accessions; s indi- 
cates Sri Lanka; u indicates Uganda; numbers above branches 
indicate bootstrap values above 50% ; letters in capital repre- 
sent genome types. 


多 基因 数据 的 合并 分 析 不 但 可 以 增加 信息 量 , 而且 能 更 可 靠 地 代表 基因 组 。 但 在 进 
行 合 并 分 析 前 ,必须 检验 不 同 基因 树 间 是 否 存在 冲突 。 通 过 PHT(partition-homogeneity test) 
( Johnson & Soltis , 1998 ) , 对 4 组 数据 的 不 同 组 合 进行 了 评估 ( 表 4) 。 结 果 表 明 , Adh 2A 
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图 3 核 Adh RAM RAWAM Æ: Adhl 树 ; A: 4dh2 树 。 种 名 后 的 数字 代表 不 同样 品 ,s 代表 斯 里 兰 卡 , u 代表 
乌干达 ; 分 支 上 的 数字 为 大 于 50% 的 自 展 支持 率 。 

Fig. 3. The maximum likelihood tree based on Adh genes. Left, Adhl; right, Adh2. Numbers after species name indicate differ- 
ent accessions; s indicates Sri Lanka; u indicates Uganda; numbers above branches indicate bootstrap values above 50% . 


表 4 不 同 数据 组 合 的 PHT 检验 结果 Adh 1 的 基因 树 间 存在 显著 的 冲突 (P = 

Table 4 Results conducted by partition-homogeneity 

test for various combinations of the data sets 0.03) , Adh2 和 其 他 片 段 的 基因 树 有 一 些 差 
异 ,但 冲突 不 显著 (P=0.11 - 0.15), Adhi, 





数据 矩阵 Data sets P 
matK vs. Adhl 1.00 =1.00)。 由 前 面 的 分 析 及 图 2、3 可 见 ， 
tin oe Adh2 和 其 他 片段 的 冲突 主要 在 O. rhi- 
eee bog zomatis 的 系统 位 置 ,为 此 我 们 对 O. rhi- 
Adhl vs. Adh2 0.03 zomatis 的 样品 进行 了 重复 扩 增 和 测序 ,系统 
a hiss 下 ie a 发 育 分 析 的 结果 仍然 不 变 , 说 明 这 一 冲突 并 
Sets Vs da Šai 非 来 自 取样 和 实验 误差 ,而 很 可 能 是 O. 
ITS vs. Adhl vs. Adh2 0.15 rhizomatis 的 Adh2 基因 片段 经 历 了 某 些 特 
MaK vs Svs Adhl vs. Adh? 0-25 殊 的 进化 事件 。 因 此 ,在 进一步 合并 DNA 


片段 数据 进行 系统 发 育 分 析 时 ,我 们 只 利用 
Adh1. matK 和 ITS 这 3 个 DNA 片段 。 


基于 合并 数据 得 到 惟一 MP 树 (图 4)(CI = 0.92, RI =0.88), 其 基本 的 拓扑 结构 和 
Adh1 单 基因 树 类 似 ,但 是 各 分 支 的 分 状 率 都 普遍 提高 ,3 个 主要 分 支 的 支持 率 均 达到 
100% ; O. officinalis-O. rhizomatis 分 支 支持 率 为 96% ,而 且 O. eichingeri 的 全 部 样品 也 以 
56% 的 支持 率 聚 合成 单独 的 分 支 ( 图 4)。 合 并 数据 的 ML 树 ( 图 4) 同样 显示 了 3 个 大 的 分 
支 ,分 别 代 表 了 不 同 的 染色 体 组 类 型 , 含 C 染色 体 组 的 全 部 二 倍 体 物 种 首先 构成 一 个 独 
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立 的 分 支 (支持 率 100% ) ,然后 与 含 B 染色 体 组 的 O. punctata 构成 姐妹 群 (支持 率 
100%), & ER EAA O. australiensis 则 与 C、B 染色 体 组 构成 的 分 支 形 成 姐妹 群 。 在 
C 染色 体 组 分 支 内 部 ,0. officinalis WMO. rhizomatis 先 组 合 在 一 起 (支持 率 9%6% ) ,再 和 0. 
eichingeri 小 分 支 构 成 姐妹 群 。 尽 管 0.，eichingeri 的 斯 里 兰 卡 样 品 和 乌 于 达 样 品 在 合并 的 
基因 树 上 构成 了 单 系 小 分 支 ,但 支持 率 较 低 (65% )。 
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ML 








MP 100 
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图 4 叶绿体 matK 基因 、nrDNA ITS 片段 和 核 基因 Adl 数据 合并 后 的 基因 树 A: 惟一 最 简约 树 ( 树 长 = 590，C1 = 
0.92, RI=0.88); A: 最 大 似 然 树 。 种 名 后 的 数字 代表 不 同样 品 ,s 代表 斯 里 兰 卡 ,u 代表 乌 干 达 ; 分 支 旁 边 数字 为 大 
于 50% 的 自 展 支 持 率 ,字母 代表 染色 体 组 类 型 。 

Fig. 4. The gene tree based on combined matK gene, ITS region and Adh1 gene data sets. Left: the single most parsimonious tree 
(tree length = 590, CI = 0.92, RI =0.88); right: the maximum likelihood tree. Numbers after species name indicate different acces- 
sions; s indicates Sri Lanka; u indicates Uganda; numbers on the branch indicate the bootstrap value above 50% ; the letters C, B, 
and E represent genome types. 


3 讨论 
3.1 染色 体 组 BC 和 下 的 系统 发 育 关系 

ERR P ,染色 体 组 B.C 和 下 的 系统 发 育 关系 一 直 存 在 争议 。Nayar (1973) 在 总 结 
稻 属 染色 体 组 间 关 系 时 指出 ,A 染色 体 组 和 C 染色 体 组 是 公认 的 基本 染色 体 组 ,B.C 和 DD 
染色 体 组 彼此 间 有 一 定 的 关系 ,E、F 染色 体 组 和 其 他 染色 体 组 的 关系 不 清楚 。Wang 等 
(1992) 采 用 核 DNA 限制 性 片段 长 度 多 态 性 (RFLP) 方 法 研究 了 稻 属 21 个 种 的 亲缘 关系 ， 
在 单独 探讨 二 倍 体 亲 缘 关 系 时 ,他 们 认为 C 染色 体 组 和 了 染色 体 组 的 关系 要 近 于 和 B 染 
色 体 组 的 关系 ,但 是 在 他 们 对 全 部 21 个 物种 构建 的 聚 类 图 上 ,E 染色 体 组 却 明 显 位 于 C 
染色 体 组 和 B 染色 体 组 共同 的 分 支 之 外 ,显示 下 染色 体 组 和 另 两 个 染色 体 组 间 较 远 的 关 
Fo Li 等 (2000) 和 凭借 基因 组 原 位 杂交 技术 (GISH) ,将 0. punctata (BB) 和 O. australiensis 
(EE) 的 总 DNA 作为 探 针 ,和 O. minuta (BBCC) 的 染色 体 组 分 别 进行 杂交 ,结果 显示 :在 前 
一 杂交 组 合 中 ,B 染色 体 组 杂交 信和 号 明显 ,C 染色 体 组 信号 微弱 ;而 在 后 一 杂交 组 合 中 ,B、 
C 染色 体 组 上 得 到 了 相同 的 杂交 信号 ,从 而 推论 B 和 C 染色 体 组 的 分 化 非常 大 ,而 下 和 
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B.C 染色 体 组 的 关系 都 很 近 ,C 和 下 染色 体 组 的 关系 要 近 于 C 和 B 染色 体 组 的 关系 。 但 
正如 Belyayev 等 (2000) 指 出 的 ,GISH 方法 在 研究 异 源 多 倍 体 的 系统 发 育 关系 上 存在 技术 
上 的 局 限 , 它 更 多 的 是 反映 重复 序列 相互 作用 的 一 般 技术 ,而 不 适用 于 探讨 系统 发 育 关 
系 。 在 我 们 目前 的 研究 中 ,3 组 DNA 片段 的 独立 和 合并 分 析 均 显示 ,C 染色 体 组 是 首先 和 
B 染色 体 组 以 100% 的 支持 率 形 成 单 系 ,而 下 染色 体 组 则 与 含 B-C 染色 体 组 的 分 支 形 成 姐 
妹 群 ,表明 B 和 CC 染色体 组 的 关系 要 近 于 它们 和 下 染 色 体 组 的 关系 (图 1 -4)。 这 个 结果 
和 Aggarwal 等 (1999) 以 及 Joshi 等 (2000) 分 别 利 用 AFLP 和 ISSR 方法 对 称 属 系统 发 育 关 系 
研究 得 出 的 结果 一 致 。 
3.2 含 C 染 色 体 组 二 倍 体 物 种 间 的 系统 发 育 关 系 

含 C 染 色 体 组 的 二 倍 体 物种 包括 分 布 在 亚洲 的 0.， officinalis #1 O. rhizomatis 以 及 间 
断 分 布 在 亚洲 和 非洲 的 O，eichingeri。 对 于 这 3 个 物种 的 界定 基本 上 没有 疑义 。0. of 
ficinalis 和 0O，rhizomatis 在 形态 上 有 明显 区 别 , O. officinalis 小 穗 较 小 ,内 释 顶 端 稍 短 ,圆锥 
花序 基部 节点 分 枝 较 多 , 等 长 , 轮 生 ,小 穗 着 生 位 置 远离 基部 的 初级 分 枝 ( Vaughan， 
1990). O. rhizomatis 则 小 穗 较 大 ,内 释 顶 端 伸 长 ,小 穗 着 生 位 置 常 靠近 圆锥 花序 的 主轴 ， 
在 圆锥 花序 基部 节点 处 缺少 轮 生 分 枝 等 (Vaughan，1990)。 另 外 由 于 杂种 不 育 , 这 两 个 种 
在 遗传 上 也 存在 一 定 差异 (Sharma & Shastry, 1965). O. eichingeri 和 以 上 2 个 种 的 主要 不 
同 在 于 没有 根 状 茎 (Vaughan，1989，1994)。 本 研究 中 ,多 基因 片段 的 独立 和 合并 分 析 均 
显示 O. officinalis WMO. rhizomatis 系统 发 育 关 系 较 近 ,而 它们 与 0，eichingeri 的 关系 则 稍 
远 (图 1-4)。 

对 于 O. eichingeri 种 内 分 化 问题 ,特别 是 斯 里 兰 卡 样品 的 分 类 地 位 一 直 不 是 非常 清 
楚 。 以 往 的 研究 多 集中 在 形态 观察 和 杂交 实验 以 及 由 此 进行 的 推断 。 例 如 ,Sharma 和 
Shastry (1965) 根据 这 个 类 群 没有 根 状 蕉 , 秆 较 细 ,叶片 窗 和 小 穗 基 部 不 着 色 等 性 状 , 将 其 
处 理 为 一 个 新 种 0. collina ; Gopalakrishnan 和 Sampath (1966) 则 认为 这 些 形态 特点 在 0. 
officinalis 变异 范围 之 内 ;通过 观察 0.collina AO. officinalis (印度) 以 及 O. eichingeri 
(乌干达 ) 分 别 杂交 时 减 数 分 裂 中 期 的 染色 体 配 对 行为 和 杂种 的 育 性 ,Sampath 和 Subra- 
manyam (1968) 认为 0. collina 是 0. officinalis WO. eichingeri 的 中 间 过 渡 类 型 ,但 倾向 将 
HEX O. eichingeri 的 亚 种 。Vaughan (1990) 在 全 面 研 究 斯 里 兰 卡 的 野生 稳 时 指出 ,分 类 
学 家 同意 O. eichingeri 是 斯 里 兰 卡 的 药 稻 复合 体 植物 的 惟一 代表 ,主要 是 受 Bor (1960) 
和 Tateoka (1962) 的 有 影响。 但 Nayar (1973) 指出 , Bor (1960) 将 O. collina 包括 在 
0 ，eichingeri 之 中 ,可 能 误 认 为 0. collina 也 是 四 倍 体 ,因为 当时 他 认为 O. eichingeri 是 四 
倍 体 。 而 Tateoka (1962) 可 能 相信 了 Bor (1960) 的 判断 ,因为 他 没有 给 出 任何 理由 就 认为 
药 称 复合 体 中 所 有 Ceylon ( 即 现在 的 斯 里 兰 卡 ) 的 材料 都 和 O. eichingeri 相同 。 同 时 ,Na- 
yar (1973) 也 提 到 , Tateoka (1962) 把 小 穗 的 尺寸 作为 区 别 O. minuta, O. officinalis , 
O. eichingeri#ll O. punctata 的 最 重要 性 状 , 那 么 按 此 标准 , 0. collina 应 该 包括 在 Tateoka 
界定 的 O. punctata 之 内 。 本 研究 结果 表明 ,尽管 斯 里 兰 卡 和 非洲 的 O，eichingeri 构成 单 
系 分 支 ,显示 彼此 较 近 的 系统 发 育 关系 ,但 它们 在 合并 的 MP 和 ML 树 上 的 支持 率 (分 别 为 
56% 和 65%) 均 显著 低 于 O. officinalis 和 0O. rhizomatis 两 个 亲缘 种 间 的 支持 率 (96%)( 图 
A) ,说 明 这 两 个 地 区 间 的 O. eichingeri 发 生 了 较 大 的 遗传 分 化 。Nayar (1973) 推测 ， 
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O. collina¥] REF O. officinalis 关系 较 近 , 它 的 形态 变异 可 能 是 由 于 周 缘 的 隔离 造成 的 ,而 
O. collina 也 可 能 是 亚洲 和 非洲 药 稳 复合 体 植物 的 过 渡 类 群 或 纽带 。 但 是 ,对 于 这 一 问 
题 的 解决 ,还 有 待 进一步 对 斯 里 兰 卡 的 称 属 植物 进行 大 量 包 括 形 态 和 地 理 在 内 的 细致 考 
察 ,从 居 群 的 角度 探讨 两 个 大 洲 间 类 群 的 遗传 关系 ,从 而 最 终 进行 可 靠 的 分 类 处 理 和 系统 
3.3 系统 发 育 研 究 中 数据 的 合并 及 其 意义 

虽然 分 子 系统 学 研究 正在 飞速 发 展 ,但 是 重建 生命 系统 发 育 树 仍然 是 具有 挑战 性 的 
难题 。 一 个 主要 的 困难 在 于 基因 的 系统 发 育 关 系 和 生物 类 群 间 的 系统 发 育 关 系 并 不 等 同 
(Wendel & Doyle, 1998; Sang, 2002)。 因 此 ,利用 单 基因 构建 的 系统 发 育 树 ,更 多 的 是 反 
映 单 基因 的 进化 历史 而 不 一 定 代 表 真 正 的 物种 系统 发 育 关系 。 为 此 ,当前 分 子 系统 发 育 
重建 的 趋势 是 通过 多 个 独立 基因 (DNA) 片 段 来 构建 多 基因 系统 发 育 树 并 进行 比较 分 析 
(Hillis, 1995; Wendel & Doyle, 1998; Ge et al. ，1999)。 因 为 多 个 相互 独立 DNA 片段 构建 
的 基因 树 更 能 反映 物种 真正 的 系统 发 育 关 系 (Miyamoto & Fitch, 1995; Sang，2002)。 而 
且 , 将 来 源 不 同 的 DNA 片段 的 序列 数据 进行 合并 处 理 可 以 有 效 地 增加 系统 发 育 信息 和 减 
少 噪 音 , 提高 系统 发 育 重建 的 准确 性 (Donoghue & Sanderson, 1992; Olmstead & Sweere, 
1994; Miyamoto & Fitch, 1995). Olmstead 和 Sweere(1994) 利 用 来 自 茄 科 17 个 种 的 2 个 叶 
绿 体 的 基因 片段 ndhF 和 rpeL 序 列 以 及 叶绿体 全 基因 组 限制 性 位 点 基因 图 谱 的 资料 , 展 
示 了 在 系统 发 育 研究 中 合并 数据 较 单个 数据 具有 更 强 的 分 辨 率 和 解释 系统 发 育 关系 的 能 
力 。 在 本 研究 中 , matK、ITS 和 Adhl 各 自 单独 的 基因 树 虽 然 不 存在 任何 冲突 ,但 是 在 分 辨 
率 和 支持 率 上 明显 不 同 。Adh1 的 基因 树 分 辨 率 最 高 ,各 分 支 平均 支持 率 为 94.9%;ITS 分 
辩 率 次 之 ,各 分 支 的 平均 支持 率 为 84% , matK 分 辩 率 最 小 ,各 分 支 平 均 支 持 率 为 76.5% 
(图 1 -3)。 当 将 这 3 组 数据 合并 后 ,不 但 分 辨 率 较 各 个 单 基因 树 要 明显 提高 ,而 且 各 分 
支 的 平均 支持 率 达 到 95.7% (图 4)。 通 过 分 析 Adhl 树 和 Adh2 树 之 间 的 冲突 ,可 以 发 现 
冲突 起 因 于 O. rhizomatiso PHT 分析 还 表明 ,如 果 将 全 部 数据 矩阵 进行 合并 将 会 掩盖 由 
Adh2 数据 造成 的 不 一 致 。 例 如 , 当 分 别 加 入 matK 或 ITS 数据 时 ,P 值 达到 0.11 或 0.15， 
均 大 于 0.05, 而 把 4 组 数据 全 部 合并 后 得 到 的 P 值 (0.25) 也 表现 出 明显 的 增 大 的 趋势 ,说 
明 hdh2 数据 中 的 噪音 在 合并 的 数据 矩阵 中 被 真实 的 系统 发 育 信 息 所 掩盖 ( 表 4) ;但 是 ， 
如 果 将 4 组 数据 强行 合并 ,得 到 的 系统 发 育 树 和 以 上 3 个 非 冲 突 数据 合并 树 比较 ,具有 冲 
突 的 分 支 的 支持 率 明显 下 降 。 因 此 ,数据 合并 之 前 进行 必要 的 一 致 性 检验 ,是 发 现 基因 树 
间 冲 突 并 决定 数据 如 何 合 并 的 关键 。 在 此 基础 上 的 合并 树 ,不 仅 提 高 了 分 辩 率 ,而 且 排 除 
了 个 别 特殊 基因 (DNA) 片 段 的 干扰 ,从 而 更 能 反映 真正 的 系统 发 育 关系 。 
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